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The recent preprint by Mazin ai. cond-mat/0212417 contains many inappropriate evaluations 
and/or criticisms on our published work [Phys. Rev. B 66, 020513 (2002) and Nature 418, 758 
(2002)]. The preprint cond-mat/0212417 vl] was submitted to Physical Review B as a comment on 
one of our papers [Phys. Rev. B 66, 020513 (2002)]. In the reviewing process, Mazin et al. have 
withdrawn many of the statements contained in cond-mat/0212417 vl, however two claims remain in 
their revised manuscript cond-mat/0212417 /3]: (1) the calculated variations of the superconducting 
energy gap within the a- or the 7r-bands are not observable in real samples due to scatterings, and 
(2) the Coulomb repulsion /i(k, k') is negligibly small between a- and 7r-states and thus should be 
approximated by a diagonal 2x2 matrix in the a and n channels. Here, we point out that the 
former does not affect the validity of our theoretical work which is for the clean limit, and that the 
latter is not correct. 



Our computational workiS for the superconducting 
properties of MgB2 is an exact implementation of the 
fully anisotropic Eliashberg formalism which was estab- 
lished more than two decades ago. The momentum de- 
pendence and the dynamical behaviors of the electron- 
phonon interactions are fully considered in this formal- 
ism. By solving the anisotropic Eliashberg equations 
correctly without any assumption on the distribution of 
the superconducting energy gap on the Fermi surface, we 
conclusively obtained the theoretical values of the super- 
conducting energy gap in MgB2 for the first time 2 .. 

Recently, Mazin et al. incorrectly stated in the origi- 
nal version of their preprint that our worki^ represents 
a computational implementation of ideas proposed by 
Liu, Mazin, and Kotus^ and that our treatment of the 
fully anisotropic Eliashberg equations integrates out all 
the phononic degrees of freedom. These statements are 
deleted in the revised version of their preprint. Mazin 
et al. also incorrectly claimed in their preprint that the 
Coulomb pseudopotential between a and it sheets might 
have been omitted erroneously in our actual computa- 
tion. This incorrect claim is also deleted in their revised 
preprint. 

The superconducting energy gap in MgB2 is shown to 
have greatly different average values on the a- and the 
7r-sheets of the Fermi surface, and it also exhibits some 
variations within the a- or the 7r-sheets£. This is the 
result of ab initio calculations of Ref . 2 carried out within 
the fully k-dependent (anisotropic) Eliashberg formalism 
in the clean limit, and no claim was made in our paper 
regarding whether the variations within the a- or the tt- 
sheets would be measurable in particular samples. We 
agree with the discussion 3,5 that these variations inside 
the a- or the 7r-sheets may be averaged out with sufficient 
strong scatterings in a sample, as also would eventually 
be the case for the difference between the a and 7r gaps. 
Since this is an impurity effect, the theoretical result in 
Ref. 2 remains valid for samples in the clean limit. We 



hope that perhaps clever experiments in the future may 
observe these variations in appropriate samples. 

Mazin et al. also claimed that the anisotropy of the 
Coulomb repulsion /i(k, k') on the Fermi surface is large 
and very important in MgB2 3 ' 5 . They further argued 
that, when viewed as a 2 x 2 matrix in the a and 7r chan- 
nels, n should be taken as a diagonal matrix (i.e., no re- 
pulsion between a and it states) as opposed to a uniform 
matrix in order to obtain physical results. This claim 
was based on their simplified estimation of /i in MgB2 by 
replacing the screened Coulomb interaction with a con- 
tact potential <5(r — r'). Their model calculation gives the 
result of /Vo- : AW ■ AW : AW = 3.1 : 2.6 : 1.4 : 1.0. 

A realistic calculation of \i requires knowing the 
screened Coulomb interaction W(r,r') which typically 
has an extent of over a bond length. Ab initio calcula- 
tions of ii can be achieved by employing the full dielec- 
tric matrij*&. This approach has been applied to MgB2 
recently and obtained a result of \i aa : /i T7r : Li a7T : 
AW = 1.75 : 2.04 : 1.61 : 1.00 following the definition 
in Ref. 5. These ab initio results show that fi is much 
less anisotropic than the model results in Ref. 5. Thus, 
the anisotropy of fi is not as important as claimed in 
Ref. 5, in particular in view of the fact that the value 
and anisotropy of A is an order of magnitude bigger in 
this system. In a simple modeling of fx, it may in fact be 
more appropriate to take it as a uniform matrix than a 
diagonal one. 

In addition, Figure 1 in Ref. 5 presents an inappro- 
priate reduction of our full k-dependent theorjiS to a 
two-band model. It is likely that the claimed discrep- 
ancies originate from their incorrect splitting of the k- 
dependent a 2 F{uj) into a 2 -Fy(w)'s for i,j = a, it. (They 
used the same frequency dependence for all components 
- cr-<7, 7T-7T, and a-iT.) A two-band model properly re- 
duced from our fully k-dependent Eliashberg formalism 
in fact reproduces well the results in Refs. 1 and 2. These 
2-band model results will be published elsewhere 8 . 
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